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Abstract 

The non-evanescent vibration, evanescent vibration and forced vibration of an oscillator attached 
to the massless spring are always discussed in general mechanics courses. In this article, we focus 
on the heavy-spring conditions. We first investigate the general situation where both viscous 
resistance and applied force are considered under the perspective of the renormalization group 
theory. Then we deal with the special case where renormalization method fails to work. We use 
analytic method to study the evanescent vibration of an oscillator attached to the heavy spring. 
Generalized orthogonal and over-complete base set is established. Two independent vibration 
modes indicate the vibrational frequency becomes bigger than simple harmonic oscillator if the 
mass of spring is considered. 
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I. INTRODUCTION 


An oscillation is a common but very important phenomenon in the physical world. If 
a physical quantity is displaced from the equilibrium a little, negative feedback may then 
lead to an oscillation. A familiar example is a simple harmonic oscillator. Also, evanescent 
vibrations and forced vibrations of an oscillator are normally focused. 1 The mass of the 
spring is neglected in models. However, the mass of the spring is unnecessarily neglected. 
In this article we try to solve the mass-spring system where the mass of the spring is not 
negligible. 

In 1979, Weinstock studied the normal modes of the oscillator motion for the oscillator 
attached to a heavy spring by virtue of the Stieltjes integral 2 . In 1994, Nunes da Silva 
obtained the normal frequencies of elastic oscillations of a particle suspended on a spring of 
non-negligible mass again under the perspective of the renormalization group theory 3 . 

A continuous spring can be regarded as a chain of many small springs coupling an equal 
amount of small masses. We then repeat mapping process by associating two consecutive 
small springs into a single one. At last, we only need to focus our attention on boundary 
effects, da Silva dealt with the problem only for the simplest situation, without the friction 
and applied force, and the advantage of this method does not emerge in the simplest case. 
In fact, we can not only find out the normal frequencies in the conservative system, but 
also obtain the specific equation of motion when external forces are acted on. In the next 
sections, we first explore the most general condition, a forced vibration with friction, with 
the help of the renormalization method. Then we deal with a special case analytically where 
the renormalization method can fail to work. We derive the equations and investigate the 
generalized orthogonality of the base set of the solution, and then obtain the normal modes 
of a damping oscillator attached to a heavy spring. 


II. THE FORCED VIBRATION WITH FRICTION CONSIDERING THE MASS 
OF THE SPRING 

Hang an uniformly distributed spring with mass m vertically. The top side is fastened 
to stable fixture and the bottom side concatenate an object with mass M as oscillator. The 
free length of spring is L and the stiffness coefficient is k. To discretize the heavy spring, 
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we view it as a series of N equal and small non-mass springs and each small spring coupling 
a concentrated object with mass m/N. The natural length and elastic constant (labeled by 
s) of each small spring are L/N and s = kN respectively. The damping on the oscillator 
M can be calculated as —bvM if the velocity of oscillator is vm and damping coefficient is b. 
Besides, we applied a time-dependent force f(t ) on the oscillator. 



FIG. 1. A discretization model for the heavy spring with natural length L, mass m and stiffness 
coefficient k when divided into N small springs coupling equal amount concentrated objects. 


The positions of objects are denoted by Xi(i — 0,1,..., N), then equations of motion can 
be set up as 


x 0 = 0, 

TTL d 2 TTL 

= s(x n _ 1 - 2x n + x n+ i) + — g (1 < n < N — 1), 
N d t z N 

(^2 

M—x n = ~s{x N - x N -i - —) + Mg - b—x N + f(t). 


( 1 ) 


We can eliminate the constant terms derived from gravity in the Eq.(l) by changing the 
coordinates appropriately. So that we can use newly defined coordinates(u*, i = 0,1, ...,N) 
to describe the motion. J.M.N. Dasilva proposed a way to approach these newly defined 
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coordinates. 3 The equations with new coordinates are written as follows, 


u 0 = 0, 
in d 2 

Ndt? Un = s( ^ Un - 1 ~ 2Un + Un+1 " > ( x < n < N ~ !), 

M—u N = -s{u N - u N _ i) - b—u N + fit), 
dt- at 


( 2 ) 


If we denote / = ^[/(t)] = ^ f f(t)e lut dt and U n = T[u n {t)\ = f u n (t)e lut dt(n = 

— OO — OO 

0,1, then the Eq.(2) can be Fourier transformed as 

Uq = 0, 

cU n = U n _ i + U n+1 (1 < n < N — 1), (3) 

CU N = U N _i + /, 
with c = 2(1 - and C = 1 + 

By Eq.(3), C 2 U2n = c(U2n—l + U2n+l) = U 2n -2 + ^U 2n + U 2n+2 ^lld cCUN = cC/jV -1 + cf ~ 
Un -2 + Un + cf. We combine two small springs into a bigger small spring. That is, with 
this mapping process, the previous (2 n — l) th and (2n) th springs now become the n th bigger 
small spring. The new Fourier transformed position function of n th spring is denoted by U n 
and U n = U 2n ■ So we have 

(c 2 - 2)U n = U n _ i + U n+ 1 (4a) 

(cC - l)Tu = Z7WT7 + cf (4b) 

Comparing Eq.(4) and Eq.(3), the equations change regularly 6 after the process of combining 

two consecutive small springs into bigger spring. So we repeat the combination to renormal- 

p-i 

If we set N = 2 P , then after p th repeat, Eq.(4b) finally becomes C^Ui = U 0 + / c^ K \ 

K=0 

P-1 


ize. 


m 


U 0 = 0 and the position function of the oscillator is M 0 sciiiator(^) = 


oo j i i c 

k ;=0 


(k) 


c(p) 


do;. 


V -1 


The next thing to do is to find out the iteration value and c^ K \ Compare between 

/y=0 

Eqs.(3) and (4) gives c^ 1 ) = c 2 — 2 and = c C — 1, so 2C (1 ) — c= c(2C — c) and more 
generally, 


p-i 


2 - c (p) = (2 C - c) JJ c ( 


k) 


Av—0 
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P -1 

Introduce 7 and let c = 2 cos 7 _ and c ^ 

K =0 

p -1 

substitution which are = 2 cos 2 p 7 and = 

K=0 


are easily accessible with this variable 


sin 2 P 7 

sin7 


. Finally, 
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c = 2 cos7 = 2(1 — Notice that > 1 , so 7 C 1 and 7 = sin7 = 


^oscillator (f) 


P -1 

co/n^ 

1 - [ -0 e ibJt duj 

2vr / C( p ) 


— OO 
OO 




2n J sin( 2 p 7 )( 2 C — c) + 2 sin 7 cos( 2 p 7 ) 

— OO 

OO 

00 J f(t)e~ lut dt 


e wt do; 


e iut du. 


2 vr J iub — Mui 2 + cot (ojyf^)uj-v Km 


(5) 


It’s the motion equation of the oscillator.' But one thing to note here is that /(f) can not 
be 0 or the solution vanishes. We deal with this condition in following part. 

OO 

In fact, if q(t) = /- f - -— -7=— -=e lu:t dcu, the motion equation of the oscil- 

v 27r -oo - iub ~ Mu + cot< ^\ff )ujdkWi ’ H 

lator can be given by the convolution of /(f) and g(t), i.e. Mosciiiator(f) = /(f) * git)- From 
the following discussion, we will know that the frequency which satisfies —iub — Mu 2 + 
cot((u = 0 is exactly the eigenvalue in no applied force condition. 


III. NO APPLIED FORCE CONDITION 

Since we can no longer use renormalization method for no applied force condition, we then 
use mathematical physics equations to study this problem. Adopt appropriate coordinates 
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as introduced before, and the problem can be analytically described as follows. 8 

u tt ~ ^ u xx = 0 {t > 0, x G [0, L)) 
u ( 0 , t) = 0 

< u x (L,t ) = -|f u tt (L,t ) - -j^u t (L, t) (6) 

u\ t= o = 0(x) 


We consider using method of separation of variables to solve this equation and we take 
u(x,t ) = X(x)exp(-iiiL as the ansatz. The equations above then become 


X"(x) + fi 2 X(x) = 0, 

X(0) = 0, 

< X'(L) = (i^+^)X(L), (7) 

X{x) = <f>(x), 

To investigate the base set of solutions, we set X p (x) and X q (x) are solutions to the equations 
above from the base set. 


x;(x) + £X r (x) = 0 ( 8 ) 

K(L) = + 4—)X t (L) (9) 

V km m 

X"(x) + ij, 2 X,(x) = 0 (10) 

X'(L) = (i M ,X= + fi 2 ,—)X,(L) (11) 

V km m 
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(9)xX q (L) — (ll)xX p (L) and we get 

'(hp — f Jj q)Xp{L)X q (L) + i(/ip — Hq) —=^=X p (L) X q (L) 


m 


\/km 


=X q {L)X'(L) - X'(L)X p {L) 


= X q (x)X'Jx) 


no 


X' p ( x ) x ' q (x)dx - X p (x)X , q (x)\ C] + / X' p (x)X' q (x)dx 

cL 


= / X,(x)d{X'(x)} - / X p (x)d{X'(x)] 


'0 


= / X q {x)X"{x) dx 


X p (x)Xq(x) &X 


( 12 ) 


= - /ip / X q (x)Xp(x)dx + /iq / Xp(x)X f/ (x)dx. 
Jo Jo 

From eq.(12) we eventually reach the following equality, 


Xq{x)X p (x)dx + 


ML 


+ i 


m (i_i p + fi q )Vkm 


X p {L)X q (L) = 0. 


(13) 


By definition, eq.(13) shows that the solutions in the base set are generalized orthogonal 5 . 
The square modulus (denoted by N 2 ) of the eigenfunctions can be calculated as 


N 2 [X p (x)}= / X 2 Jx)dx + 


ML 


m 


+ i- 


2/j,pVkm 


X 2 P (L) 


Xp(x)=sinfi p x LI cos/ipLsm/ipL ML 9 

-sm2/ipL H-1-sm a p L 

° " 2 m 


(14) 


Eq.(9) 


L 

2 


ML 

2m 


4/L 

sin 2 n p L. 


2/ip 


So the solution to the Eq.(6) can be written as u(x,t) = zl n sin/i n a;exp(— i/i n L-y^t) and 

the expansion coefficient can be determined from the initial conditions. 0(x) and ^(x) 
could be expanded as 4>(x) = -P n sin/i n x, //?(x) = —i /i n Ly^|Q n sin/i n a; where 

0(L)sin/i n L 


Pn 


Qr, 


Jq (j)(x)sin/.i n xdx + ^ + i 


m 2 ii n y/km 


N 2 [X n (x )] 

Jo* iiL\/fV'(^)sin/i n ^ + 


— + i 

Ln 2 fi n y/km 




(15) 


iV 2 [X n (x)] 

Let = a?„ + /5Q„, we then have 

0(x) = + /3Q n )smn n x, 


ip{x) = —i VV n LW — (aP n + /3Q n )smn r 

V m 


(16) 


x. 
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Compare Eq.(16) with 0(» = Yl n P n sm/.i n x, ip(x) = -i J2 n Vn p \ £-Q n sm/jL n x, 


P n [a(l - i/JL n L\ /—) - 1] + Qn[/3(1 - in n L\,—) + ifx n L\ —] = 0, 

V m V m V m 

P n [a(l + i{i n L\ /—) - 1] + Q n \fi{ 1 + in n L\ —) - i/J. n L\.I—] = 0. 

V m V m \ m 


( 17 ) 


Hence, a = 


A n — 


1 


and [5 = 




1 ii/Zn Z/ 
r L 


N^[X n (x)](l ± i^L^) LA 


. Finally can be given as 


(x) =F -0(^)]sin / n ri a;dx + + i- ^-=][(j)(L) =p ^(L^sin/^L 

m 2 \x n \Jkm 


(18) 


This non-unique expansion also suggests that the generalized orthogonality of base set 
{X n {;x)} is over-complete. 

The eigenvalue equation reveals as cot /iL = i^== + fx^A (eigenvalue // won’t be 0). 
Notice that —JJ is also eigenvalue if // is eigenvalue. The corresponding expansion coefficients 


have the relation P(—/x) = —P(/x) and Q(—/x) = —Q(/u) due to the Eq.(15). Considering 


sin(— /xx) = —sin/xx, 


y, P n siri/i n .T = y [p n siri/i n x + P n siri/i n 2 )] = y [P n sin/z n x + c.c.], 

n Re(/i n )>0 Re(/i n )>0 

- z y fx n L\ — Q n sh\fi n x = -i[ y [jx n L\ — Q n sm/j, n x] - c.c.]. 

V m V m 

n Re(n n )>0 


(19) 


Eq.(19) verihes that both cj)(x) and ip( x ) are pure real when ansatz u(x, t ) = AiSin/i n xexp(— in n L\/ A 

f~k 

is taken. Of course u(x,t ) = ) [A n sm/j, n xexp(—in n L\ — t ) +c.c.] is also pure real. We 

V m 

Re(£tn)>0 

can obtain two independent eigen-vibration modes from this result. Let /i n = n 9 , then 

F 


smu n xexp(—iu n L\ —t) 

V m 

= sin(£ n - iC n )icexp[-i(f n - i( n )L\ — t] 

V 777/ 

= (sin^xcosh^a: — icos£ n £sinh(( n a;)exp(— C n L\ — t)[cos(£ n L\ —t) — isin(£ n L\ —0] 

V m V m V rn 

= [sin&jZCOshCnacos (£ n L\ —t) - cos£ n xsmh( n xsm(€ n L\ —t)]exp(-( n L,J—t)~ 

\ m V m V m 

i[cos£ n xsmh( n xcos(£ n L\ — t) + sin£ n xcosh£n£sin(£ n LA/— t)]exp(-( n LW — t). 

V m V m V m 


( 20 ) 













e ^ 10 cos(0 e ‘t 5 cos(0 




FIG. 2. Let t be dimensionless with r = y^/£ n L, the eigen-vibration mode of the oscillator can 
be plotted under two cases with Cn/£n = 0.1 and Cn/£n = 0.2. The oscillator vibrates with damped 
amplitude. 

The two independent vibration modes are given by 

mode l:[sin£ n xcosh() n xcos(£ n Z/'\/— t) — cos£ n xsmh( n xsm(£ n L\ — t)]exp(—( n L\ —t), 

V m V m, V m 

mode 2:\cos£ n xsmhQ l xcos(£ n L\ /— t) + sm£ n xcosh( n xsm(£ n L\ — t)]exp(— ( n L\ — t ). 

V m V m V m 

The oscillator (x = L) vibrates with damped amplitude (as is shown in Fig.2) in both modes, 
which is reasonable. 

For zero-friction case (b = 0), the solutions to eigenvalue equation cot fiL = i-r== + fi^- 
are pure real, i.e. ( = 0. Then two independent vibration modes become 

k 


model: sin£ n xcos(£ n L\/— t ), 

*' m 

k 


mode 2: sin^xsin^Li/— t ). 

*' m 


With cot £ n L = we can also obtain > \J JJ- So the vibrational frequency 

becomes bigger than simple harmonic oscillator if the mass of spring is considered. 10 

Summarize the result, the solution to Eq.(6) is u(x,t) = Yl n AiSin/r n xexp(— ifi n L^J^t) 
where eigenvalue fi n is given by cot fi n L = = + p-nq/p- Expansion coefficient A n and 

modulus square are given by Eqs.(18) and (14). 


IV. CONCLUSION 


In this article, we detailedly studied the vibration of spring oscillator when the mass of 
the spring can’t be neglect. Damped oscillation and forced vibration are especially focused. 
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For general condition, oscillation with friction and applied force, renormalization method is 


employed to obtain the equation of the motion. Renormalization method shows superiority 


when there is applied force f(t) exerts on the oscillator. We also investigate the damping 
vibration without applied force with theory of partial differential equations. For this spe¬ 
cific boundary condition, the generalized orthogonality of base set is studied. The motion 
equation of this condition is also given. We discussed the characters of the eigenvalue and 
the expansion coefficient and the discussion verifies the validity of the solution. 
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k f f{t)e~ iult dt 

— oo 


k f dt 

— OO 


e i(—cj)i d(—ca), 


oo oo 


oo 


oo 



pure real. 

Take infinitesimal spring with length dx, stiffness coefficient ^ and mass ^-m. Suppose its 
amount of deformation is d u. Then its equation of motion in infinitesimal time dt can be given 


10 







by m ■ 4tt according to Newton’s second law. 


9 It’s hard to find all /j ,s analytically from cot fiL = £=== + but we can use perturbation 
to investigate since imaginary part in RHS of eigenvalue equation is feeble in generally 
underdamping conditions. The zero-order approximation, ^°)tan//°) L = for zero friction 
situation (6 = 0), is discussed in previous work. 9 Obviously, solutions are real numbers and 
come in pairs (—/^t 0 ) is also solution if satisfies zero-order eigenvalue equation). Perturba¬ 
tion calculation shows the first order modification is negative pure imaginary number and the 
paired zero-order eigenvalues it// 0 ) share the same first order modification. It coincides with 

the previous analysis that —/Z is eigenvalue if n is eigenvalue. 
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If we neglect the mass of the spring, the frequency £ n L\/ — degrades, lim £ n L J — = \ tj, 


m—>0 


which can be proved by using lim cotx — ^ = 0. 

x—>0 x 
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